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ABSTRACT Complete Middle East respiratory syndrome coronavirus (MERS-CoV) genome sequences were obtained from nasal 
swabs of dromedary camels sampled in the Kingdom of Saudi Arabia through direct analysis of nucleic acid extracts or following 
virus isolation in cell culture. Consensus dromedary MERS-CoV genome sequences were the same with either template source 
and identical to published human MERS-CoV sequences. However, in contrast to individual human cases, where only clonal 
genomic sequences are reported, detailed population analyses revealed the presence of more than one genomic variant in indi- 
vidual dromedaries. If humans are truly infected only with clonal virus populations, we must entertain a model for interspecies 
transmission of MERS-CoV wherein only specific genotypes are capable of passing bottleneck selection. 

IMPORTANCE In most cases of Middle East respiratory syndrome (MERS), the route for human infection with the causative agent, 
MERS coronavirus (MERS-CoV), is unknown. Antibodies to and viral nucleic acids of MERS-CoV have been found in dromedar- 
ies, suggesting the possibility that they may serve as a reservoir or vector for human infection. However, neither whole viral 
genomic sequence nor infectious virus has been isolated from dromedaries or other animals in Saudi Arabia. Here, we report 
recovery of MERS-CoV from nasal swabs of dromedaries, demonstrate that MERS-CoV whole-genome consensus sequences 
from dromedaries and humans are indistinguishable, and show that dromedaries can be simultaneously infected vsdth more than 
one MERS-CoV. Together with data indicating widespread dromedary infection in the Kingdom of Saudi Arabia, these findings 
support the plausibility of a role for dromedaries in human infection. 
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Two hundred twelve cases of Middle East respiratory syndrome 
(MERS), 88 of them fatal, have been reported since April 2012 
(1). Although examples of human-to-human transmission have 
been identified, the origin of infection with the causative agent, 
MERS coronavirus (MERS-CoV), is unexplained in the majority 
of cases (2). Serologic evidence of infection in dromedary camels 
(DC) and, more recently, the detection of viral nucleic acid in DC, 
particularly in juvenile DC, suggest the possibility that DC may 
serve as a reservoir or vector for human infection (3-12). How- 
ever, there are as yet no published analyses of complete MERS- 
CoV genomic sequences or virus isolation from DC in the King- 
dom of Saudi Arabia (KSA). 

In a collaborative effort between the Center for Infection and 
Immunity in the Mailman School of Public Health at Columbia 
University and the Mammals Research Chair, Department of Zo- 



ology, College of Science, King Saud University, a mobile labora- 
tory was established in Saudi Arabia to investigate the possible role 
of DC, other domestic animals, and wildlife in the transmission of 
MERS-CoV through molecular and serological analyses. In a pre- 
vious publication, we reported detection of high loads of MERS- 
CoV nucleic acid in nasal swabs from DC (10). Here, we describe 
MERS-CoV complete genome sequencing, detailed phylogenetic 
analyses, and the recovery of live virus through culture. 

Reverse transcription-PCR (RT-PCR) assays of nasal swab 
samples demonstrated the presence of MERS-CoV RNA in DC at 
a high prevalence in KSA (10). Sequence analysis of products rep- 
resenting three regions of the MERS-CoV genome revealed iden- 
tity over approximately 3,000 nucleotides (nt) with human 
MERS-CoV sequences. To determine whether this identity ex- 
tended across larger regions of the MERS-CoV genome, we pur- 
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TABLE 1 High-throughput sequencing of MERS-CoV from dromedary camels in Saudi Arabia 



Sequencing approach 





Ion Torrent 






lUumina 








Random-primed 


Preamplified RT-PCR 


Random-primed 


Random-primed 


Preamplified RT-PCR 


Random-primed 




cDNA (swab 


fragments (swab 


cDNA (cukure 


cDNA (swab 


fragments (swab 


cDNA (culture 


Sample' 


extract)" 


extract)' 


extract)'^ 


extract)" 


extract)*" 


extract)' 


KSA-363n 


Yes 


Yes 


Yes 


Yes 


Yes 


NA 


KSA-378n 


Yes 


NA"* 


Yes 


Yes 


NA 


NA 


KSA-376n 


Yes 


Yes 


NA 


Yes 


Yes 


NA 


KSA-344r 


Yes 


Yes 


NA 


Yes 


Yes 


NA 


KSA-409n 


Yes 


Yes 


NA 


Yes 


Yes 


NA 



" Sequencing library prepared from random-primed cDNA that was generated from total nucleic acid extract of nasal swab sample. 

^ Sequencing library prepared from pooled overlapping RT-PCR fragments (14) amplified from total nucleic acid extract of nasal swab sample. 

' Sequencing library prepared from random-primed cDNA that was generated from total nucleic acid extract of cell culture supernatant from infected Vero cells. 

^ NA, not applicable. 

' n, nasal swab sample; r, rectal swab sample. 



sued whole-genome sequencing using the Ion Torrent and lUu- 
mina platforms employing as the template random-primed cDNA 
libraries and pools of PGR products based on primers that repre- 
sented published human MERS-CoV genomic sequence. 

Raw Ion Torrent and lUumina data from 5 DC were assembled 
against MERS-CoV scaffolds available from GenBank. No 
platform-dependent differences were apparent; thus, sequence 
data were combined and used to assemble consensus sequences 
for each sample. The specific processing of individual samples is 
summarized in Table 1. Consensus full-genome sequences of 
MERS-CoV from DC were obtained for samples KSA-363-Taif- 
21, KSA-378-Taif-36, and KSA-376-Taif-34 (10). Partial genomes 
were obtained for samples KSA-344-Taif-2 and KSA-409-Tabuk- 
26. 

Two additional full genomic sequences were generated entirely 
by overlapping direct RT-PCR amplification of random-primed 
cDNA generated from total nucleic acid extract of Arabian DC 
nasal swab samples KSA-503-Taif-45 and KSA-505-Taif-47. 
These samples were collected in January 2014 in Taif from a 
1 -year-old imported African and a 1 -year-old Arabian breed of 
DC, respectively. 

Complete genomic sequence was also obtained for virus cul- 
tured from two DC nasal swab samples. Vero E6 cells were inoc- 
ulated with sterile filtered nasal swab/viral transport medium 
(VTM) samples (KSA-363-Taif-21 and KSA-378-Taif-36) or a 
rectal swab/VTM sample (KSA-371-Taif-29 [10]). Viral prolifer- 
ation was monitored by real-time "upstream-of-E" (UpE) PCR 
after 48 h and 66 h in comparison to residual inoculum measured 
after removal of the inoculum and washing of the cells at 2.5 h 
postinfection {t = 0) (Fig. lA). Virus growth was observed with 
the two nasal swab samples but not with the rectal swab sample. 
Total nucleic acid extracts obtained from the 48 -h samples were 
subjected to random sequencing on the Ion Torrent platform, 
yielding fiiU-length genomic sequence. No differences were ob- 
served in the consensus sequences obtained using template from 
extracts of nasal swabs or cultured virus. 

The five consensus fiiU-genome sequences (KSA-363, -378, 
-376, -503, and -505) were aligned to other genome-length human 
MERS-CoV sequences available in GenBank. Analysis of the five 
consensus sequences confirmed earlier work with short PCR 
products obtained from DC (10) that suggested that DC in Saudi 
Arabia harbor the same virus that causes MERS in humans 
(Fig. IB). 



However, detailed inspection of the multiple sequence align- 
ment indicated that our DC sequence assemblies showed frequent 
lAUC codes for two-base wobbles in positions where high diver- 
gence between human MERS-CoV sequences is observed (Fig. IC, 
i). To ensure that the appearance of multiple sequence variants in 
individual DC samples was genuine rather than a sequencing ar- 
tifact, we amplified representative regions by specific PCR, cloned 
the products, and sequenced individual clones by the dideoxy 
chain termination method. Alignment of the clone-derived se- 
quences confirmed the presence of multiple sequence species in 
several individual DC samples (Fig. IC, ii and iii). Sequence diver- 
sity in sample KSA-363 decreased over a period of 48 h in culture 
(see Fig. SI in the supplemental material). 

The role of DC in human MERS-CoV infection is unclear. 
Studies throughout the Middle East conducted independently by 
several research teams have described antibodies to MERS-CoV in 
DC (3-12). In recent work in KSA, we found antibodies to MERS- 
CoV in 95% of adult DC and MERS-CoV sequences in 35% of 
juvenile DC (10). However, direct exposure to DC is only rarely 
reported in human cases. Furthermore, there are no published 
reports of MERS-CoV virus isolation from DC and only a single 
near-full-genome sequence from an African DC that indicates 
that viruses related to human MERS-CoV circulate in African DC 
(12). Here, we confirm that DC may harbor infectious virus and 
that whole-genome consensus sequences obtained from nasal iso- 
lates align with whole-genome sequences recovered from humans. 

Our analysis of whole-MERS-CoV-genome sequences recov- 
ered from DC revealed the presence of sequence variants within 
single samples (also known as quasispecies [13]). One amino acid 
in the spike protein (A520S, corresponding to nucleotides [nt] 
23013 to 23015 in GenBank accession no. JX869059) was changed 
within the receptor-binding domain; however, all other changes 
occurred outside the receptor-binding domain. It is unclear 
whether this has any functional implications. In other viral sys- 
tems, genetic diversity has been linked to pathogenicity and 
shown to enable adaptation to new environments such as those 
associated with movement into new hosts (13, 14). No sequence 
variants have been described in individual human MERS-CoV 
samples. Whether this means that only consensus sequences are 
reported or that human sequences truly represent clonal virus 
populations within individual cases cannot be discerned from 
published data. If the latter, we must entertain a model for inter- 
species transmission, wherein only specific genotypes, which may 



2 mBio' mbio.astm.org 



May/June 2014 Volumes Issues e01146-14 



Dromedary MERS-CoV Genomic Sequence 



,0^ 



B 




■ KSA 363 
i KSA 378 



0 hrs 



48 hrs 



66 hrs 



|— AI-Hasa-1 7 KSA 201 3(KF600647) 
-AI-Hasa-21 KSA2013(KF600634) 
- AI-Hasa-1 6 KSA 201 3 (KF6C0644) 
AI-Hasa-2 KSA201 3 (KF1 86566) 
■AI-Hasa-152013(KF600645) 
AI-Hasa-1 KSA2013(KF186567) 

Riyadh-9 KSA 201 3 (KJ1 56869) 

England 2 201 3 
L Riyadh-2KSA2012(KF600652) 
Buraidah-1 KSA2013(KF600630) 
KSA-CAIVIEL-376 



881— 



— France UAE2013(KF745068) 

- Hafr-AI-Batin-2 KSA2013 (KJ1 5691 0) 

Riyadh-14KSA2013(KJ156934) 

J- Hafr-AI-Batin-1 2013(KF600628) 
1 Madinah-3 KSA 201 3 (KJ 1 56950;kJ 1 5691 6) 

rQatar3 2013(KF961221) 
lool Qatar42013(KF961222) 
^Riyadh-1 KSA2012(KF600612) 
lool Bisha-1 KSA 2012 (KF600620) 

KSA-CAMEL-503 

KSA-CAMEL-363 

KSA-CAMEL-50S 



- KSA-CAMEL-378 

■Munich UAB 201 3 {KF1 92507) 
England Qatar 201 2 (KC667074) 
- Riyadh^ KSA201 3 (KJ1 56952) 

- Riyadh-3 KSA 201 3 (KF60061 3) 



64L 



- Riyadh-5KSA2013(KJ156944) 
-Wadi-Ad-Dawaslr-1KSA2013(KJ156881) 
-Taif-1 KSA2013(KJ156949) 



J EMC Bisha 201 2 (JX369059) 

n Jordan Al-Zarqa 201 2 (KC7761 74) 

NRCE-HKU205-Egypt{KJ477102) 

FIG 1 (A) Real-time PGR analysis of cell culture supernatant after inoculation of Vero cells with nasal swab samples KSA-363 and KSA-378. (B) Phylogenetic 
analysis of MERS-CoV sequences from dromedary camels in Saudi Arabia and other genome-length MERS-CoV sequences available on 7 AprU 2014. GenBank 
accession numbers are given in parentheses for each sequence {England2 sequence is available at http://www.hpa.org.uk/Topics/InfectiousDiseases/Infections 
AZ/MERSGoV/respPartialgeneticsequenceofnovelcoronavirus/); bootstrap values of >60o/o indicate statistical support for the respective nodes; the scale bar 
indicates the number of substitutions/site. (G) (i) Glippings from the multiple MERS-GoV sequence alignment indicating sequence variation among human 
MERS-GoV sequences and potential variation of sequences within individual DG samples (indicated by x). (ii) Sequences obtained by direct sequencing of PGR 
products from the same region, (iii) Sequence analysis of individual clones generated from the PGR amplification products. Two PGRs were performed, including 
nt 24190 to 24300 and 24510 to 2530. (Gontinued on following page) 



not be present in every infected DC, are capable of passing bottle- 
neck selection. Such a model would not abrogate a role for host 
susceptibility in infection and disease but might provide insights 
into the rarity of human cases of MERS. 

Nucleic acid extraction, high-throughput sequencing, and 
PGR. Total nucleic acids from nasal swab, rectal swab or cell cul- 
ture supernatant samples were extracted on a QiaCube with Cador 
reagent kits (Qiagen, Hilden, Germany). Superscript III and ran- 
dom hexamer primers were used to generate cDNA preparations 
(Life Technologies, Carlsbad, CA, USA). Second-strand cDNA 
synthesis for high-throughput sequencing was carried out by ran- 
dom primer extension with Klenow enzyme (New England Bio- 
labs, Ipswich, MA, USA). 

High-throughput sequencing was performed in parallel on 
random-primed cDNA preparations and MERS-CoV-enriched 
PCR product pools. To enrich for MERS-CoV sequences from 
total nucleic acid extracts, PCR amplifications employing a set of 
overlapping PCR primers spanning the whole genome in approx- 
imately 2.0- to 2.5-kb fragments were performed as described 



(15). PCR products were pooled and sequenced on both the Ion 
Torrent and lUumina platforms. 

Sequencing on the lUumina HiSeq 2500 platform (lUumina, 
San Diego, CA, USA) resulted in an average of approximately 20 to 
50 million reads per sample. cDNA preparations were sheared 
(E210 sonicator; Covaris, Woburn, MA, USA) for an average frag- 
ment size of 200 bp and added to Agencourt AMPure XP beads 
(Beckman Coulter, Brea, CA, USA) for purification, and libraries 
were prepared with Kapa high-throughput library preparation 
kits (Kapa Biosystems, Wilmington, MA, USA). Sequencing was 
performed using a read length of 100 nt, followed by an indepen- 
dent read of the 6-nt bar code. Samples were demultiplexed using 
Illumina-supplied CASAVA software and exported as FastQ files. 
More than 90% of lUumina reads passed the Q30 filter. Demulti- 
plexed FastQ files were mapped against GenBank scaffolds 
(KF600620 and KF186567) with Bowtie 2 mapper 2.0.6 (http: 
//bowtie-bio. sourceforge.net [16]). Sequencing on the Ion Tor- 
rent PGM platform was performed with Ion PGM Sequencing 200 
kits on Ion 318 chips (Life Technologies), yielding on average 1.5 
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FIG 1 (Continued) 



to 2.5 million reads per sample with a mean length of approxi- 
mately 166 nt. cDNA preparations were sheared (Ion Shear Plus 
kit; Life Technologies) for an average fragment size of 200 bp and 
added to Agencourt AMPure XP beads (Beckman Coulter, Brea, 
CA, USA) for purification, libraries were prepared with Kapa li- 
brary preparation/Ion Torrent series kits (Kapa), and emulsion 
PGR was performed with Ion PGM Template OT2 200 kits (Life 
Technologies). Ion Torrent reads were demultiplexed and ex- 
ported as FastQ files by the Ion Torrent PGM software. After bar 
code and adaptor trimming, length filtering, masking of low- 
complexity regions, and subtraction of ribosomal and host se- 
quences, reads were mapped as described for lUumina data. Con- 
sensus sequences from mapping assemblies were generated by 
using SAMtools/BGFtools 0.1.19 software (http://samtools 
.sourceforge.net [17]). 

Based on available sequence information, a set of 25 nested 
consensus primer sets were designed to generate overlapping PGR 
products of approximately 1.3 kb that comprise the full genome 
(see Table SI in the supplemental material). Random-primed 
cDNA was PGR amplified with individual primer pairs and Am- 
pliTaq Gold (Life Technologies). The PGR products were purified 
by agarose gel electrophoresis and QIAquick gel extraction kits 
(Qiagen) and subsequently sequenced on both strands by the 
dideoxynucleotide chain termination method (GeneWiz, South 



Plainfield, NJ, USA). Products from selected PGRs were also 
cloned into pGEM-T Easy plasmid vector (Life Technologies), 
and 16 individual clones were dideoxy sequenced in order to as- 
sess clonal sequence diversity. 

Quantitative real-time PGR used OneStep Real-Time qPGR 
buffer (Life Technologies) and UpE primer/probes (18). 

Virus isolation. One hundred fifty microliters of nasal swab in 
universal virus transport medium (Becton, Dickinson, Franklin 
Lakes, NJ, USA) was filtered (0.45 ijum; Millipore, Billerica, MA, 
USA), and the filtrate was inoculated on Vero E6 cells grown to 
semiconfluence in T25 culture flasks with Dulbecco modified Ea- 
gle medium (DMEM)-10% fetal calf serum. The inoculum was 
removed after 2.5 h, cells were gently rinsed, and fresh medium 
was added (T-0) . Supernatant was tested for MERS-GoV by quan- 
titative real-time PGR after 48 h (T-48), and supernatant as well as 
cell homogenate harvested at 66 h postinfection (T-66). 

Bioinformatics and phylogenetic analysis. Sequence data 
were analyzed using software packages Geneious (Biomatters, 
Auckland, New Zealand), MEGA (http://www.megasoftware.net 
[19]), and Wisconsin GGG (Accelrys Inc., San Diego, GA). Phy- 
logenetic analysis was performed by the neighbor-joining method 
implemented in MEGA 5.2, running 1,000 pseudoreplicate anal- 
yses to assess statistical support. 
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Nucleotide sequence accession numbers. Full genomic se- 
quences of MERS-CoV from DC were deposited in GenBank un- 
der the indicated accession numbers: samples KSA-363-Taif-21, 
KJ713298; KSA-378-Taif-36, KJ713296; KSA-376-Taif-34, 
KJ713299; KSA-503-Taif-45, KJ713297; and KSA-505-Taif-47, 
KJ7 13295. 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi:10.1128/mBio.01146-14/-/DCSupplemental. 

Figure Si, PDF file, 0.5 MB. 

Table Si, DOCX file, 0.1 MB. 
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